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ABSTRACT 

Young stars form in molecular cores, which are dense condensations within 
molecular clouds. We have searched for molecular cores traced by 13 CO J = 1 — > 
emission in the Taurus molecular cloud and studied their properties. Our data 
set has a spatial dynamic range (the ratio of linear map size to the pixel size) of 
about 1000 and spectrally resolved velocity information, which together allow a 
systematic examination of the distribution and dynamic state of 13 CO cores in a 
large contiguous region. We use empirical fit to the CO and C0 2 ice to correct 
for depletion of gas-phase CO. The 13 CO core mass function ( 13 CO CMF) can be 
fitted better with a log-normal function than with a power law function. We also 
extract cores and calculate the 13 CO CMF based on the integrated intensity of 
13 CO and the CMF from 2MASS. We demonstrate that there exists core blending, 
i.e. combined structures that are incoherent in velocity but continuous in column 
density. 

The core velocity dispersion (CVD), which is the variance of the core velocity 
difference 5v, exhibits a power-law behavior as a function of the apparent sepa- 
ration L: CVD (km/s) oc L(pc) ' 7 . This is similar to Larson's law for the velocity 
dispersion of the gas. The peak velocities of 13 CO cores do not deviate from the 
centroid velocities of the ambient 12 CO gas by more than half of the line width. 
The low velocity dispersion among cores, the close similarity between CVD and 
Larson's law, and the small separation between core centroid velocities and the 
ambient gas all suggest that molecular cores condense out of the diffuse gas with- 
out additional energy from star formation or significant impact from converging 
flows. 
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Introduction 



Most young stars are found in dense molecular cores (IMcKee &: Ostrikerll2007l ). There is 
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One of the primary outcomes of these studies is the core mass function (CMF). Because of 
the still unknown physical origin of the stellar initial mass function (IMF) and its signif- 
icanc e, emphasis has b een placed on the possible connection between the CMF and the 
IMF jAlves et al.ll2007h . 



The construction of an analytic form of the CMF from observational data has largely 
focused on two functional forms, power law and log-normal. The majority of pas t studies 
claim to find a power law CMF, the shape of which resembles the Salpeter IMF fjSalpeter 
19551 ) (The upper mass limit of the original Salpeter IMF is 10M Q ) 
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dlogM 



oc M" 7 , 7 = 1.35, -0.4 < log(M/M Q ) < 1.0 . 



Su ch a power law CM F was found in millimeter continuum maps of the p Ophiuchus region 
by iMotte et al.l ( 119981 ) with the IRAM 30- m eter telescope and a serie s of subsequent studies 



in the Serpens regi on ( Testi fc Sargentl 1998), the p Op hiuchus region (Johnstone et al. 2000 
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Johnstone fc Ballvl 120061 ) and the Pipe nebula (lAlves et al.l 120071 ). 
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Reid fc Wilson! ( 12006bl ) studied the CMF in 11 star- forming regions and find an average 
power law index of 7 = 1.4 ± 0.1. Recently, dust emission observations of the Aquila 
rift complex with Herschel [j] re veal a power law mass function with 7 = 1.5 ± 0.2, for 
M > 2 M Q (jKonyves et al.ll2010l ). which is also consistent with the Salpeter IMF. A similar 
power law index is found in some molecular emission studi es. For example, the CMF obtained 



from a C 18 study in the S140 region has 7 = 1.1 ± 0.2 ( Ikeda fc Kitamurall2011af) . On the 



other hand, there are also studies that find a flatter CMF. iKramer et al.l (119981 ) studied seven 
molecular clouds L1457, MCLD 126.6+24.5, NGC 1499 SW, Orion B South, S140, M 1 7 SW 
and NGC 7538 in 13 CO and C 18 0, and find 7 to be between 0.6 and 0.8. iLi et all J2007h 



x http: / /www.esa.int /SPECIALS /Herschel/index.html 
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studied cores in the Orion molecular cloud traced by sub-millimeter continuum and found a 
power law with 7 = 0.15 ± 0.21. These findings of s mall 7 are in the minori ty and do not 
seem to be a special result of spectroscopic mapping. iRathborne et al.l (120091 ) obtained the 
CMF from an extinction map and used complementary C 18 observations to examine the 
effects of blending of cores in dust maps. They claimed to find a CMF with 7 similar to that 
of the Salpeter IMF. 

At a first glance, a similarity between the CMF and Salpeter IMF suggests a constant 
star formation efficiency, which is independent of the core mass. It is crucial to note, however, 
th ese studies are examini ng structures on vastly different scales of size, mass, and density. 
In iReid fc Wilson! (j2006al ). for example, the mass of the cores ranges from about 0.1 M & to 
1.6 x 10 4 Mq. Due to the large distances of many targeted regions, any "cores" over about 
500 M & are certainly unresolved, with m any showing signs of much evo lved star formation , 
such as water masers (jWang et al.l 120061 ) and/or compact HII regions (IHofner et al.ll2002l ). 
The observed similarity between CMF and Salpeter IMF may be explained equally well by 
self-similar cloud structures as well as a constant star formation efficiency. 



Some observations (e.g., Enoch et al.l ( 120081 )) suggest a log- normal form for the CMF 



(in the mass range of 0.1 M Q < M < 10 M 
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This is reminiscent of the IChabrierl (120031 ) IMF, which is also of log-normal form. Theoreti- 
cally, if the core mass depends on n quantities which are random variables, th e CMF would be 
log-n ormal when n is large, i.e., the core formation processes are complicated (lAdams &: Fatuzzo 
19961 ). This is a result of the central limit theorem . A log- normal distribution also arises 
naturally from isothermal turbulence ( )Larsonlll973l ). 



It is thus of great interest to distin guish the two form s of th e CMF and obtain the key 
parameters associated with each form. ISwift fc Beaumontl (120101 ) show that a large sample 
with many cores is needed to differentiate these two forms. Furthermore, we also emphasize 
here the critical need to obtain a large sample of cores in spectroscopic data. Overlapping 
cores along the same line of sight can only be separated using resolved velocity information. 
It is important to evaluate the effect of such accidental alignment on the derived CMF. A 
Nyquist sampled continuous spectroscopic map is also essential for studying the dynamical 
characteristics of star forming regions, such as the Core Velocity Dispersion (CVD = (dv 2 ) 1 ! 2 
see section 4.3), of the whole core sample in one star forming region. 



The Taurus molecular cloud is a nearby (with a distance of 140 pc. iTorres et al.l (120091 )) 
low-mass star-forming region. In this work, we obtain a sample of cores in the 13 CO data 
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cube of this region. We study the properties of 13 CO cores in detail and compare them with 
those found in the dust extinction map of the same region. We first briefly describe the data 
in §21 we present the methods used to find 13 CO cores in §31 we present the observed 13 CO 
CMF and CVD in §U we discuss the implications of our observations in §0 In the final 
section we present our conclusions. 



The Data 



In this work, we use 12 CO and 13 CO data in the form of a (x, y,v) cube of the Taurus 
molecular cloud as observed with the 13.7 m FCRA O telescope (INaravanan et al.l 120081 ) 



and the 2MASS extinction map of the same region (jPineda et al.l 120101 ). The 12 CO and 
13 CO lines were observed simultaneously between 2003 and 2005. The map is centered at 
a(2000.0) = 04 /l 32 m 44.6 s , 5(2000.0) = 24°25'13.08", with an area of ~ 98 deg 2 . The FWHM 
beam width of the telescope is 45" at 11 5 GHz. The angular s pacing (pixel size) of the 
resampled on the fly (OTF) data is 20" ([Goldsmith et al.l 120081 ) . which corresponds to a 
physical scale of ~ 0.014 pc at a distance of D — 140 pc. There are 80 and 76 velocity 
channels in the 12 CO and 13 CO data cube, respectively. The width of a velocity channel is 
Kh = 0.266 km/s. The extinction map has a pixel size of about 5 times that of the CO data 
cube with, of course, no velocity information. 

The probability density function (PDF) of 13 CO data and that of 2MASS extinction 
data are shown in Figure [TJ from which the difference between the two data set can be seen. 
The normalization of 13 CO integrated intensity and that of 2MASS extinction are chosen so 
that they roughly correspond to the same column density of molecular hydrogen. 



3. Core Extraction 



Cores have been empirically defined as regions with concentrated, enhanced intensity 
in a data cube or a map. We empirically assume an ellipsoidal shape for a core and use the 
FINDCLUMPS tool in the CUPID package, which is a part of the star link software^. We 
have tried two methods in the FINDCLUMPS tool, GAUSSCLUMPS Jstutzki fc Guesten 



1990|) and CLUMPFIND f lWilliams et al.lll994f ). GAUSSCLUMPS searches for an ellipsoid 
with Gaussian density profile around the brightest peak and subsequently subtracts it from 
the data. It then continues the process with the core-removed data, iterating successively 



2 http://starlink.jach. hawaii.edu/starlink/ 
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until a terminating criterion is reached. CLUMPFIND identifies cores by drawing enclos- 
ing contours around intensity peaks without assuming the shape of cores a priori. Unlike 
GAUSSCLUMPS, CLUMPFIND cannot deconvolve overlapping cores. Our subsequent anal- 
ysis is based on GAUSSCLUMPS and we discuss the caveats of CLUMPFIND at the end of 
this section. 

For fitted ellipsoids of revolution, the core radius is defined as the geometrical mean of 
the semi-major and the semi-minor axes 

R = (-Rmax-Rmin) ^ • (3) 

We take the observed size R as a typical scale of a core. 

Following the instructions for CUPID, we first subtract the background by using FIND- 
BACK. Some studies of the Taurus molecular cloud find a characteristic length scale of 
about 0.5pc ~ 2pc, at which self-similarity breaks down and gravity becomes important 



( IBlitz fc Williams! 119971 ). This length scale corresponds roughly to 35 ~ 140 pixels in the 
13 CO data cube, and to 7 ~ 28 pixels in the extinction map. We set the smoothing scale to 
127 pixels in each axis for the 13 CO data cube, and 25 pixels for the extinction map, both 
are close to the upper value of the scale at which the self-similarity breaks down. 

After background subtraction, we use the GAUSSCLUMPS method of the FINDCLUMPS 
tool to fit Gaussian components in the 13 CO data cube , which are identified as 13 CO cores 



if their peak intensity is larger than a threshold (see e.g. ICurtis Sz Richer! ( 120101 )). which is 
set to five times of the rms noise (parameter RMS). Since the data cube is large, running 
the 3D-Gaussian fitting on the whole data cube is time-consuming. Furthermore, the back- 
ground differs substantially between different parts of the Taurus cloud. The noise level also 
differs among different regions. The data cube is divided into several sections in the x-y 
plane as shown in figure [2] to make the sizes of data sets suitable to handle, and to minimize 
the variation of the background and the noise level. The fitting parameters are given in 
tabled] In practice, two adjacent regions are arranged to have an overlapping stripe of 100 
pixels in width to avoid the cores touching boundaries being missed. In a final step, the 
resultant catalog is checked and the duplicate cores are removed. 

FINDCLUMPS outputs the total intensity, T tot , through the summation of all pixels in 
each fitted core. We then calculate the mass of a core based on T tot . The central frequency 
of the 13 CO J = 1 — > line v is 110.2 GHz. The column density of 13 CO in the upper-level 



(J = 1) can be expressed as (IWilson et al.ll2009l ) 
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where k is Boltzmann's constant, h is Planck's constant, c is the speed of light, A u \ is the 
spontaneous decay rate from the upper level to the lower level, and Tj, is the brightness 
temperature. A convenient form of this equation is 
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Th e tota l 13 CO column density N tot is related to the upper level column density N u through 
(O l2002h 

Ntot, 13 CO = fufrfbN u l3 CO . (6) 

In the equation above, the level correction factor f u can be calculated analytically under the 
assumption of local thermal equilibrium (LTE) as 
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where g u is the statistical weight of the upper-level. T cx is the excitation temperature 
and Q(T P ^) = k T P ^/B P is the LTE partition function, where B e is the rotational constant 
(jTennysonl 120051 ) . A convenient form of the partition function is then Q(T CX ) m T CX /2.76K. 



The correction factor for opacity is defined as 

/ n 3 dv 
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where ri3 is the opacity of the 13 CO transition and Tb g is the background temperature, 
assumed to be 2.7K. 

The 13 CO opacity is estimated as follows. Assuming equal excitation temperatures for 
the two isotopologues, the ratio of the brightness temperature of 12 CO to that of 13 CO can 
be written as 

T b , 12 1 - 
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Assuming t\i 3> 1, the opacity of CO can be written as 
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The excitation temperature T ex is obtained from the 12 CO intensity. First, the maximum 
intensity in the spectrum of each pixel is found. This quantity is denoted by T max . The 
excitation temperature is calculated by solving the following equation 
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where h, k and v are Planck's constant, Boltzmann's constant, and the central frequency of 
i2co J = 1 line (115.27 GHz), respectively. 

The total number of H2 molecules in a core as traced by 13 CO is 

N tot is co (DA) 2 



J H 2 



ri3 



CO]/[H 2 ] 



(13) 



where the distance of the Taurus molecular cloud, D = 140 pc, the pixel size of the data cube, 
A = 20", the 13 CO to H 2 abundance ratio [ 13 CO]/[H 2 ] is taken to be 1.7x 10" 6 JFrerking et al. 



19821 ). Calculation of the mass of the 13 CO-traced core is then straightforward 

M = /3m H2 SH 2 , 



(14) 



where /3 = 1.39 c onverts the hydrogen mass to total mass taking into account of the abun- 
dance of helium (IWilson Sz Roodlll994h . 



We then correct for the depletion of 13 CO by using the recipe of IWhittet et al.1 (120071 ). 
who measured the column density of CO and C0 2 ices toward a sample of stars located 
behind the Taurus molecular cloud and obtained the following empirical relationship 
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The column density of the depleted CO is given by iV(CO)* c ° e tal = iV(CO) icc + iV(C0 2 ) icc , 
the total column density of CO is then iV(CO) total = iV(CO)t tal + N(CO ) ms . Using the 
abundance ratio of 13 CO to 12 CO, fisco /xgcQ = 1/70 Jstahl fc WilsorJ [l992h. the depletion 
correction for Taurus can be calculated. iPineda et al.l f )2010l ) obtained a good linear corre- 
lation between the extinction and the total 13 CO column density (gas plus ice) using such a 
correction. We first calculate the average extinction in the 13 CO core region, which corre- 
sponds to a column density of 13 CO and 13 C0 2 ice. Since a 13 CO core occupies only a part 
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of the velocity channels, the ratio of the 13 CO integrated intensity and total intensity of the 
13 CO core region, / cor e/-f total is calculated and the mass correction for a 13 CO core is 

AM = fi3 CO/12co nR 2 N(CO)\T l ■ (17) 

where R is the typical scale of a 13 C0 core defined in equation [3j This recipe achieves good 
depletion correction at the arc-minute scale. 

The FWHM line width AVfwhm is calculated from the velocity dispersion of each 13 CO 
core, Av, which is the standard deviation of the velocity value about centroid velocity, 
weighted by the corresponding pixel data value. 



AVfwhm = 2^2 In 2 Aw . (It 



The virial mass is also estimated from Av ( iBertoldi fc McKedll992l ) 



5Av 2 R 

vir = ~~ ' ( ^ 

This equation describes a balance between the self gravity and combined thermal and non- 
thermal motions of a 13 CO core, neglecting external pressure and the magnetic field. A core 
is considered gravitationally bound when its mass exceeds M v - a . 

In dealing wi th the extinction m ap, the hydrogen column density is estimated from 



optical extinction ( iGuver &: Ozelj 120091 ) 



= 2.2 x 10 21 . (20) 

cm 1 \ mag / 

We have applied GAUSSCLUMPS with the fitting parameters listed in table [TJ to the 
Taurus region to extract cores. We have also applied CLUMPFIND with the fitting parame- 
ters listed in table |2j In both cases, a background was first subtracted with the FINDBACK 
procedure using the same parameters listed in the tables. The cores fitted by GAUSS- 
CLUMPS are shown in Figure [2] and those by CLUMPFIND in Figure [31 There are many 
more smaller cores found by CLUMPFIND, especially in the relatively diffuse region. This is 
attributable to CLUMPFIND 's requirement to assign one pixel to one particular core with- 
out 'splitting' a pixel in possible overlapping configurations. We rely on the cores found by 
GAUSSCLUMPS in the subsequent analysis, since CLUMPFI ND cannot sp l it ove rlapping 



cores properly. In a study of core extraction with CLUMPFIND. iPineda et al.l (120091 ) suggest 
not to blindly use CLUMPFIND to derive mass functions in any "crowded" situation. 

Parameters of 3D 13 CO cores and those cores found in the smoothed 13 CO data cube 
can be found in tables 0] and [5j respectively. In table [61 we give parameters of the cores 
derived from the 2MASS extinction map. 
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Results 



We have performed a test on GAUSSCLUMPS by fitting to simulated data. We put 
100 Gaussian components with random sizes and random locations into an empty simulation 
cube (an empty 3D array) with the size of region 1 in Figure[2j White noise (with a ~ 0.1 K, 
a typical value across the Taurus region.) is also added to the whole cube as a background. 
This simulated data cube is then fitted by using the GAUSSCLUMPS tool. As a result, we 
find that for cores with peak intensity higher than 0.7K (7a), the distribution of the cores 
found by GAUSSCLUMPS is similar to the simulated distribution. Figu re [4] gives numerical 
resul ts. For lower masses, GAUSSCLUMPS tends to find more cores (IStutzki fc Guesten 
19901 ). Similar simulations have been done for the 2D data ( 13 CO total intensity data and 
extinction data), where a threshold of 5a is found to ensure agreement between measured 
and simulated distributions. 

From the GAUSSCLUMPS fitting, we select cores with peak intensity higher than 7a ~ 
0.7 K for 13 CO data cube based on the simulation mentioned above, 5a ~ 1.3 K-km/s for 
13 CO total intensity, and 5a ~ 1.5 mag for the extinction map, where a is the variance of 
data (the RMS). Other fitting parameters can be found in tables [H El and[3j 

Some filamentary structures (ellipses with a large axis ratio) were found in both the 
extinction map and the 13 CO total intensity map. They comprise only a small fraction 
of the total mass (less than 10%). They do not show up in the fitting to the 13 CO data 
cube, which means they are not "coherent", meaning that, the velocity variation within 
the structure cannot be described by a Ga ussian in velocity space. Since the stability of 
a filament is different from that of a core (ILombardi fc Bertml l200ll ). we filter out those 
gaussian components that have an large axis ratio (major/minor >10) in the analysis. We 



thus obtain a sample of 765 13 C O cores. This is t 



cores for a contiguous region (cf. iTatematsu et al. 
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2009 



20001 : llkeda et all 120071 
2011bh . 



Figure [2] shows the 13 CO integrated intensity map overlayed with the cores obtained by 
Gaussian fitting to the 13 CO data cube. Most of the regions with 13 CO emission are found 
to contain cores. However, despite the strong 13 CO emission in region 11 and 5, no 13 CO 
core is found in the former and only one 13 CO core is found in the latter. This is due to 
the rapid variation of intensity between velocity channels, which is clear in figures and El 
which are the channel maps of region 11 and region 5, res pectively. Gas in regio n 11 is likely 
to be affected by the young protostellar cluster within it (jSchmalzl et al.ll2010l ). In contrast 
to these two regions, the 13 CO emission in all other regions containing 13 CO cores is found 
to change gradually in velocity (see, e.g. figure [7]). 
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Fig. 1. — The PDFs of 13 CO integrated intensity (dot-dashed line) and 2MASS extinc- 
tion (solid line) are shown. The horizontal axis shows the 13 CO integrated intensity and 
the 2MASS extinction (normalized to values that correspond to similar column density of 
molecular hydrogen), while the vertical axis shows the number of pixels with different values 
(left: 13 CO integrated intensity; right: 2MASS extinction). 
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Fig. 2. — Image showing the total intensity of the CO emission in Taurus region. The 
data are divided into 18 regions for 13 CO core fitting. The cores found by fitting gaussians to 
the 13 CO data cube (x, y, v) using GAUSSCLUMPS are shown as green ellipsoids. Despite 
the strong 13 CO emission, no 13 CO cores are found in region 11, and only one 13 CO core is 
found in region 5. This can be understood by looking at the channel maps (see discussion 
in section 3). 




Fig. 3.— The cores found with CLUMPFIND routine are overlayed on the 13 CO total 
intensity map. It is clear from this figure that the cores are "crowded". In this case, 
CLUMPFIND cannot extract cores properly. 



T(K) 



Fig. 4. — Core extraction from simulated data for 13 CO by using GAUSSCLUMPS. The 
dashed line shows the distribution as a function of temperature of the simulated cores, while 
the solid line describes the cores extracted. When the core peak intensity is higher than 
7 a ~ 0.7 K, the difference is small. For peak intensity lower than 7a, GAUSSCLUMPS 
tends to extract more cores than are actually present. In this work only 13 CO cores with 
peak intensity higher than 7a are included. 
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Fig. 5. — 13 CO channel maps of the region 11 in figure [2j The intensity changes drastically 
between channels, which is in stark contrast with the region shown in figure [7] and suggests 
velocity incoherence in region 11. This is the reason why no 13 CO core can be fitted by a 
three dimensional (x,y,v) gaussian in this region. 
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Fig. 6. — Channel maps of region 5 in figure El The intensity here also changes drastically 
between channels. Only one 13 CO core is found in this region. 



- 16 - 




Fig. 7. — Channel maps of region 10 in figure [2] The main features are coherent through 
multiple channels and in consequence many 13 CO cores are found in this region. 
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The distributions of the peak optical depth and the mean H2 density of the 13 C0 cores 
are shown in Figure |HJ The mean H 2 density is calculated by dividing the total number of 
H 2 in a 13 C0 core by its volume, assuming its size along the line of sight is the same as its 
typical size (geometrical mean of major and minor, 2 times the typical radius in equation [3]). 
Most of the 13 CO cores have a typical J = 1 — > optical depth of ~ 0.8 and a typical mean 
H2 density of ~ 2000 cm -3 . The recipes that we have adopted for optical depth correction 
and excitation state correction are adequate for cores in Taurus. 
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Fig. 8. — The distribution of the peak optical depth (left panel) and mean H 2 density (right 
panel) of the 13 CO cores. 

By looking at the motion of the 13 CO cores, we find a group of 13 CO cores with system- 
atically different line of sight velocities (see the lower corner of figure IT6|) . which suggests 
these 13 CO cores may not belong to the main cloud. After excluding these, we are left with 
588 13 CO cores. 

For comparison, we also use GAUSSCLUMPS to fit cores in the 13 CO integrated inten- 
sity map with the procedures used for core fitting in the 13 CO data cube. Figure [9] shows 
the 13 CO integrated intensity map overlaid with cores obtained by Gaussian fit of the 13 CO 
integrated intensity map, which is a 2D core fitting. Compared with the 3D core fitting to 
the 13 CO data cube, there are cores found in the 2D core fitting within the regions having 
relatively abrupt, significant velocity variations seen in the 13 CO data. This indicates that 
velocity is essential for reasonably defining cores in a complex cloud such as Taurus. The 
physical plausibility of requiring that a core have a relatively well-defined velocity in addi- 
tional to a compact spatial structure eliminates many spurious cores that are identified in 
the 2D data alone. 

Figure [10] shows dust extinction map overlayed with cores obtained by Gaussian fit to 
the 2MASS extinction map. The 200" spatial resolution of 2MASS extinction map is about 
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5 times coarser than that of 13 CO map. The extinction cores are generally larger than those 
found in 13 CO data. We perform an experiment by smoothing the 13 CO data to the same 
spatial resolution as 2MASS extinction and perform the core fitting. The fitted cores are 
also generally larger than those found in the original 13 CO data, confirming that the cores 
found in the 2MASS data are generally unresolved even at the modest distance of Taurus 
(see table [5]). 



4.1. Mass distribution of the cores 

We analyze the CMF in two functional forms, power law (equation [T]) and log-normal 
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Therefore, we adopt a Monte Carlo approach similar to that used in iLi et al. 
first generate a random sample of cores with one of the above distributions. For a power 
law distribution, there are three parameters, M max , M min , and 7, while for a log-normal 
distribution, there are two, a and \x. The cumulative distribution C(Mj), % — 1, ...,n of the 
random sample is then fitted to that of the sample found in the data Co(Mj), % — 1, n by 
minimizing a \ 2 function defined by 



X 



[C(M t )-C (M t )]- 



(21) 



The resulting mass function of the 13 CO cores can be seen in figure [TTJ This mass 
distribution can be fitted much better by a log-normal function than by a power law, which 
only fits the higher mass end. The fitting results are similar for cores based on 2D fitting. 
The mass distribution that is based on 2D fitting of the 13 CO integrated intensity can be 
fitted better with a lognormal function. The mass function of 2MASS extinction cores can 
also be fitted better with a log-normal function than with a power law (figure [T2|) . Thus 
log-normal is a better representation of the mass distribution of cores than power law in the 
case of cores represented by 13 CO, 13 CO total intensity, and 2MASS extinction. 

To examine the reliability and the completene ss of the cores found, we calculate the 
minimum detectable mass for a core of certain size (ILi et al.ll2007l ) 



M = M, 



point 



x VN x VM 



(22) 



where M po i nt is the minimum detectable mass of a point source, and N and M are the number 
of pixels and number of velocity channels occupied by a core, respectively. We found this 
minimum detection mass is about 0.03 M Q for a core as large as the smoothing scale. There 
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is also another way to estimate the completeness (IPineda et al.ll2009l ). with which similar 
result is obtained. All the cores found in the 13 CO data cube are larger than this minimum 
mass. For the extinction map, the minimum detectable mass of a core is found to be about 
0.8 Mq. 



We used bootstrap (see e.g. iPress et al.lll992l ) to estimate the uncertainty of the fitting. 
The resulting uncertainty is extremely small when the sample size is large. This means 
that the statistical uncertainty in our Monte-Carlo type fitting procedures is negligible. The 
difference between the fitted function and the data has to be systematic. For this reason, 
we do not show the numerical values of the uncertainties. 



4.2. Energy State of "CO Cores 

We study the energy state of 13 CO cores by analyzing the mass and line width together, 
which can only be accomplished through spectral maps. The properties of the 10 most 
massive 13 CO cores are listed in table H] for the 13 CO data cube. 



We investigate the mass-radius relation of the 13 CO cores, which can be seen in figures 
and [T41 The c ore stability is studied by calculating the Bonnor-Ebert mass-radius relation 



(ILi et al.ll2007l ). in which the Bonnor-Ebert mass is obtained from integration of the density 



profile of the Bonnor-Ebert sphere 

P(0 = 1 + (g/2.25)°-s ' (23) 
where the dimensionless radius £ is defined as 

i = r^^Gp c jvl , (24) 



where p c is the central density and v s = y/ kT eq / jjuiu is the effective sound speed. The 
equivalent temperature T pq is related to the FWHM line width AVfwhm with following 



equation (ILi et al.l 120071 ) 



_ m H (3AV F 2 WHM ) 
lcq ~ 81n(2)lfc " (25j 

In our sample, the FWHM line width of the 13 CO cores ranges from 0.5 km/s to 1.7 km/s. 
When £ is larger than £ max = 6.5, there is no longer any stable solution. With £ max and the 
size of a 13 CO core, the central density p c is determined. The critical mass is then calculated 
by integration to the density profile (equation [23!) . Since the upper limit and the lower limit 
of the integration are constants, it is straightforward to show that the Bonnor-Ebert mass 
M BE oc r. 
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A 13 CO core having mass greater than its Bonnor-Ebert mass would be hydrostatically 
unstable. As can be seen in figure 12, no such 13 CO cores (assuming a FWHM line width of 
1 km/s) are found. T his is in contra st to cores found in Orion, which are mostly supercritical 
(and thus unstable) (ILi et al.ll2007l ). 



We have studied the mass distribution of nearly gravitationally bound 13 CO cores (with 
virial mass to mass ratio M vir /M < 2) and find that they can be well fitted with a log-normal 
function with /i = 2.35, a = 0.52 (see figure [T5"j) . The fact that a log-normal function fits 
the nearly gravitationally bound cores better than an Salpeter IMF type power law suggests 
that the mass conversion efficiency is NOT constant for all 13 CO cores. 

We also examine the virial mass of the 13 CO cores according to Equation [19j Among 
the 13 CO cores, there are 83 cores out of 588 that are nearly gravitationally bound (with 
the virial parameter M vir /M < 2, see equation ITS]) . The pressure from the surrounding 
materials also ser ves as confinement in ad dition to self-gravity, reducing the virial mass of 



a 13 CO core (e.g. iKainulainen et al.l (120111 ) ) . We can define an effective mass corresponding 



to the pressure by calculating the square of the effective velocity dispersion 

AV c 2 ff = — = fcrext Pext , (26) 

Pcore Pcore 

where T ext is the temperature of the ambient medium. An effective mass is then defined 
analogously to the virial mass 

Mes = fAV&R = fkT^R Pext = 4 g (f_\ ( R \ (T±*\ (PcxA M& ^ (27) 

G m H2 G p cor e ' V 5 / \0-lpc) V 10K / VPcore/ & 

where / is a dimensionless factor. For cores of the same mass, the pressure is a more 
important factor for the confinement of the less dense cores. The external pressure may 
be negligible for the confinement of a dense core, for example, considering a core with 
Pcxt/Pcorc ~ 1/100, the effective mass would be about 0.05 M , which is negligible for a 
core having a total mass > 1 M Q . The pressure is then unimportant for the confinement of 
massive cores, but is important for a core having total mass < 0.1 M . 



4.3. The Motion of Cores 



Previo us studies suggest that CO traces t he bulk motion of dense gas as well as NH 3 

I with the 



Benson et al. 


1998; 


Kirk et al. 


2010) 



fitted centroid radial velocity coded in color. In general, the separation between the velocity 
of the core and that of the ambient gas is small. As seen in figure [T71 the core centroid 
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velocities are also mostly similar, except for cores in the lower right corner. This region 
differs systematically from the main Taurus cloud and may not be physically connected to 
it. After excluding these cores, there are 588 cores (out of 765 cores) left. We plot the 
velocity difference 5v vs. the apparent separation L of these cores in figure [T8"l 

Between the 588 cores there are 588 x (588 — l)/2 unique pairs, which are shown as 
the background density distribution in figure [18j The velocity difference 5v shows a bimodal 
distribution. We then calculate the dispersion of the velocity difference and plot the core 
velocity dispersion (CVD = (5V 2 ) 1 / 2 ). The relation between the CVD and the apparent sepa- 
ration between cores can be fitted with a power law, yielding CVD (km/s) = 0.2L(pc) a7 + 0.2 
for L between and 10 pc, while those points having L > 10 pc show large scatter, likely 
caused by under-sampling due to the finite size of the cloud. The mean value of the CVD for 
separation L > 10 pc is 1.18 km/s. The bimodal distribution in the velocity difference can 
be understood as defining two spatial scales. One corresponds to the typical size of a "core 
cluster" , about 4 pc, while the other one shows the typical core cluster separation, about 8 
pc. 

To evaluate the significance of the bimodal distribution, we constructed a simple two- 
core-cluster model with 500 cores, with each cluster having a radius R = 7 pc. The clusters 
are separated by 9 pc. A uniform spatial distribution of cores is assumed for both clusters. 
The velocity of a core is defined by a Gaussian with variance a, given by (a / cr max ) = (Dj i?) a5 , 
where D is the distance of the core from the center of the cluster it belongs to, and <r max is 
the maximum a at the edge of a cluster, taken here to be 1.3 km/s. In practice, the line of 
sight velocity of a core is taken to be a random number from a gaussian distribution with 
variance There is an assumed systematic difference between the line of sight velocities 

of the two clusters, which is 1.8 km/s in our model. We are able to reproduce the major 
features of the observed CVD plot (figure [TBI) with this simple model. Rigorously speaking, 
the velocity of cores in each cluster do not follow the exact form of Larson's law, because 
there is a center. However, the effect of such a distinction is negligible. We have tried to 
generate a sample of cores in which the velocity of a core is still gaussian distributed whose 
variance depends on its distance to the previously generated core, denoted as /, with the 
same relation a oc I ' 5 , and get similar results to those shown in figure [JBJ The two length 
scales, i.e. the size of a core cluster, and the separation between the clusters are clearly seen 
in figure [19j The relation between the CVD and the apparent separation in the L = ~ 10 
pc region can be fitted by a power law of the form CVD (km/s) = 0.13L(pc) 0,7 + 0.2. In 
L > 10 pc region, CVD can be fitted by a linear function of slope 0.04. The mean value of 
the core velocity difference in this region is 0.93 km/s. We conclude that the main features 
of the observed CVD vs L relationship and the observed 5v are reproduced reasonably well 
by this simple demonstration. 
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Fig. 9. — The cores found through a 2D Gaussian fitting of the 13 CO total intensity map 
are overlaid on the 13 CO total intensity map of the whole Taurus region. The red dashed 
lines are the coordinate lines of constant right ascension and declination. 



-23 - 



mag 



28 



26 



° 24 



22 



20 



04 h 56 m 



16 




48 



40 



32 



24 



16 



oU 



EA 



Fig. 10. — The cores found through a 2D Gaussian fitting of the extinction map are overlaid 
on the dust extinction map of the whole Taurus region. They are generally larger than the 
13 CO cores. The red dashed lines are the coordinate lines of constant right ascension and 
declination. 
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Fig. 11. — The mass distribution of the 13 C0 cores found by GAUSSCLUMP fitting to the 
13 CO data cube. The left panel shows the best fit power law distribution, while the right 
panel shows the best fit log-normal distribution, with the latter giving a better fit over the 
full mass range. The minimum detection mass for the 13 CO cube is 0.03 M Q , indicated by 
the vertical dash-dotted line. We have fitted the mass function to 13 CO cores more massive 
than the minimum detection mass. The fitting results are : (a) power-law distribution with 
dN/dlogM oc M -0 - 53 , M max = 40.03 M , M min = 0.74 M . (b) log-normal distribution 
with fi = 0.95, a = 1.15. 
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Fig. 12.— The mass distribution of the cores found by GAUSSCLUMP fitting to the 13 CO 
total intensity map (upper row) and to the extinction map (lower row). The panels on the left 
are fitted with power law distributions, while those on the right are fitted with log-normal 
distributions. The minimum detection masses are 0.03 M & and 0.8 M Q for the 13 CO data 
cube and the extinction map, respectively, and are indicated by vertical dash-dotted lines. 
The fitted results are: (a) dN/dlogM oc M" a3 °, M max = 106.48 M , M min = 0.13 M . 
(b) log-normal distribution with \i = -0.62, a = 0.70. (c) dN/ d log M oc M" ai , M max = 
92.0 M Q , M min = 1.2 M . (d) log-normal distribution with // = 2.04, a = 0.77. The cores 
from 13 CO total intensity and the extinction cores are both fitted better with a log-normal 
distribution than with a power law. 
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Fig. 13.— The mass-size relation of the 13 CO cores found by GAUSSCLUMPS in the 
13 CO data cube. The Bonnor-Ebert mass is indicated by the dashed lines. The upper line 
corresponds to FWHM linewidth of 1.5 km/s and the lower line corresponds to 1.0 km/s.The 
mass-size relation is best fitted asMoc R 2 6 indicated by the dash-dot line. 13 CO cores with 
virial mass to mass ratio M vir /M < 2 are denoted by blue circles; those with 2 < M vir /M < 5 
are denoted by green diamonds; other 13 CO cores with M vir /M > 5 are denoted by black 
diamonds. 
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Fig. 14. — The relation between the mass and size of the 13 CO cores found in the extinction 
map (a) and in the 13 CO total intensity map (b) with GAUSSCLUMPS. The mass-size 
relationship of extinction cores is best fitted with a power law M cx R 2 ' , while 13 CO total 
intensity cores can be fitted with a power law M cx R 2A . The Bonnor-Ebert mass is indicated 
by the dashed lines. The upper line corresponds to FWHM linewidth of 1.5 km/s and the 
lower line corresponds to 1.0 km/s. 
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Fig. 15. — The mass distribution of those nearly gravitationally bound 13 C0 cores with 
virial mass to mass ratio M v ; r /M < 2. The fitting results are (a) dN/dlogM oc M a20 , 
M max = 40.04 M Q , M min = 3.57 M Q . (b) log-normal distribution with fi = 2.35, a = 0.52. 
The log-normal distribution gives a good fit over the entire mass range that is sampled. 
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Fig. 16. — Overlay of the cores identified using the 3D 13 CO data cube on the 12 CO total 
intensity map with the centroid velocity coded as shown in the color bar on right. In most 
parts of Taurus, the core centroid velocities are similar. However, at the lower right corner, 
the centroid velocities of the cores differ systematically from those in other regions, which 
may suggest a different origin or location of these cores. In the study of the core velocity 
dispersion, we exclude these cores. 



-30 - 




Fig. 17. — Overlay of the cores found in the 3D 13 CO data cube on the 12 CO total intensity 
map. The difference between centroid velocity of the cores and the local 12 CO peak velocity 
is indicated by the color (coded in the color bar on right) and is generally small (< 1 km/s). 
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Fig. 18. — Plot of the core velocity difference, 5v vs. the apparent separation L of cores and 
the core velocity dispersion (CVD= (Sv 2 ) 1 ^ 2 ). The background is a image of the number 
of data points in the 5v — L plane, with grey scale (density bar at top of figure) showing 
the density of points. Contours of the density distribution in this plot are also shown. The 
green diamonds represent the variance of the velocity difference in each separation bin. The 
background density distribution shows two group of points, which are also evident by looking 
at the contours. These two groups define two length scales. One corresponds to the typical 
size of a "core cluster", about 4 pc, while the other one shows the "core cluster" separation 
to be about 8 pc. 



-32 - 



Discussion 



The shape of the core mass function has been for some time under scrutiny due to its 
possible relevance to the origin of stellar initial mass function. We found that a log-normal 
distribution best represents the 13 CO CMF in Taurus, through fitting the 3D 13 CO data 
cube. 

We also searched for and studied the cores identified in the 2MASS extinction map. 
For extinction cores, a power law function is a better fit than a log-normal function. It 
is important to note the significant differences between 3D fitting and 2D fitting. In the 
former case, velocity information is invoked so that overlapping cores can be separated. We 
have compared the fitting to the 3D data cube and that to the 2D total intensity map 
(figure [201) . There are some low-density cores that do not emerge in the fitting to the 3D 
data cube. There are overlapping cores, which can be separated in the 3D fitting by using 
the velocity information (an example of the spectrum is shown in figure [2"Tj) . This is a direct 
demonstration of the peril of solely relying on total column density maps (such as dust 
continuum) to obtain core properties. In the light of these considerations, we believe that 
the CMF derived from extinction maps needs to be treated with some caution for interpreting 
the origin of the IMF, even though it may have a similar power-law form as that of the IMF. 

The large sample size and substantial spatial dynamic range of the Taurus 13 CO core 
sample allow us to reconstruct the core velocity dispersion (CVD) for the first time. The 
CVD exhibits a power-law behavior as a function of the apparent separation L between 
cores for L < 10 pc (see figure [T5j) . This is similar to Larson's law for the velocity dispersion 
of the gas. The peak velocities of 13 CO cores do not deviate from the centroid velocities 
of ambient 12 CO by more than half the 12 CO line wi dth. The small y elocity differences 
between dense and diffuse gas have also been noted by iKirk et al.l (120071 ) for Perseus cores. 
Simulation of core formation under the influence of converging flows sug gest that massive 
cores exhibit relatively small line widths compared to less massive ones (IGong &: Ostriker 
20111 ). We do not see this trend in our work, i.e. there is no apparent correlation between the 
mass and the l ine w idth or between the mass and the temperature (see figures [2"2l and [2"3"j) . 
Padoan et al.l (120011 ) found smaller line width in denser gas, which is not seen in Taurus 
sample either. These results suggest that dense cores condense out of the more diffuse gas 
without additional energy input from sources, such as protostars or converging flows. 



In recent simulations of core formation 
formation 



le-g. 



(e.g. 



Gong fc Ostrikerl 1201 ll ) or star cluster 



Offner et al.l 120091 ). there is sufficient information to produce a CVD plot of 



the simulated core samples. Since CVD is sensitive to the dynamic history of core formation, 
we encourage the simulators to perform such analysis to facilitate direct comparison between 
theoretical models and observations. 
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6. Conclusions 

We have studied the 13 CO cores identified within the Taurus molecular cloud using a 
100 degree 2 13 CO 7=1—7-0 map of this region. The spatial resolution of 0.014 pc and 
the velocity resolution of 0.266 km/s facilitate a detailed study of the physical conditions of 
13 CO cores in Taurus. The spatial dynamic range (the ratio of linear map size to the Nyquist 
sampling interval) of 1000 of our data set allows examination of the collective motions of 
13 CO cores and their relationship to their surroundings. We have found that the velocity 
information helps to exclude cores which we consider to be spurious. 

Our conclusions regarding the extraction of 13 CO cores and their properties are: 

1) Velocity information is essential in resolving overlapping cores, allowing better reso- 
lution of cores and a more accurate determination of the CMF. 

2) The mass function of the 3D(x,y,v) 13 CO cores can be fitted better with a log- normal 
function (fi = 0.95 and a = 1.15) than with a power law function. For cores represented by 
13 CO total intensity and extinction, a log- normal distribution is also a better representation 
of the mass distribution than is a power law. There is no simple relation between the Taurus 
13 CO CMF and the stellar IMF. 

3) No 13 CO cores are found to have mass greater than the critical Bonnor-Ebert mass, 
in contrast to cores in Orion. 

4) Only 10% of 13 CO cores are approximately bound by their own gravity (with the 
virial parameter M vir /M < 2). These can be well fitted with a log-normal mass function. 
External pressure plausibly plays a significant role in confining the 13 CO cores with small 
density contrast to the surrounding medium. 

5) In Taurus, the relation between core velocity dispersion (CVD= (Sv 2 ) 1 ^ 2 ) and the 
apparent separation between cores L can be fitted with a power law of the form CVD (km/s) 
= 0.2L(pc) 7 + 0.2 in the < L < 10 pc region, similar to the Larson's law, with a median 
value of 0.78 km/s. 

6) The observed CVD is reproduced by using a simple two- core- cluster model, in which 
there are two core clusters with radius 7 pc and a separation 9 pc between these two clusters. 

7) The low velocity dispersion among cores, the close similarity between CVD and 
Larson's law, and the small difference between core centroid velocities and the ambient 
diffuse gas all suggest that dense cores condense out of the diffuse gas without additional 
energy input and are consistent with an ISM evolution picture without significant feedback 
from star formation or significant impact from converging flows. 
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8) The CVD can be an important diagnostic of the core dynamics and the cloud evo- 
lution. We encourage simulators to provide comparable information based on their calcula- 
tions. 
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Fig. 19. — Results from numerically generated sample taken from two spherical clusters of 
cores, each having a 7 pc radius and with a separation of 9 pc between the centers of two 
clusters. The same statistical treatment as that applied to the actual data (figure IT8l) has 
been employed. The two length scales, i.e. the size of a core cluster, and the separation of 
the clusters are clearly seen. The points with apparent separation > L < 10 pc can be 
fitted with a power law of CVD (km/s)= 0.13L(pc) 0,7 + 0.2. The horizontal line shows the 
mean 5v value, 0.93 km/s, of points with L > 10 pc. 
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Fig. 20. — Fitting to a patch within region 10. The spectrum at the location of the red circle 
in the upper right of panel (a) is shown in figure [2TJ (a) The total intensity map of 13 CO 
cores found in the data cube of region 10 (3D).(b) 13 CO cores found in the total intensity of 
the original data of region 10 (2D). More cores are found in the 2D fitting, especially in the 
relatively diffuse regions. These cores contain large internal velocity variations. 
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Fig. 21. — An example spectrum containing contributions from two CO cores, which are 
clearly separated in velocity. The location where this spectrum was observed is indicated by 
the red circle in panel (a) of figure |20j 
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Fig. 22. — The mass-line width relation of the cores found by GAUSSCLUMPS. There is 
no apparent correlation between the mass and the line width. 
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Fig. 23. — The mass-temperature relation of the cores found by GAUSSCLUMPS. There is 
no apparent correlation between these two quantities. 
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Table 1: Parameters used in the GAUSSCLUMPS fitting of the 13 C0 data cube and 13 C0 
total intensity map. 



r ai ametei 


Value 


vv vv id i n 


z 


WMIN 


0.01 


MAXSKIP 


100 


THRESH 


5 


NPAD 


100 


MAXBAD 


0.05 


VELORES 


2 


MODELLIM 


0.05 


MINPIX 


16 


FWHMBEAM 


2 


MAXCLUMPS 


2147483647 


MAXNF 


200 



WWIDTH is the ratio of the width of the weighting function ( which is a Gaussian function ) to the width 
of the initial guessed Gaussian function. 

WMIN specifies the minimum weight. Pixels with weight smaller than this value are not included in the 
fitting process. 

MAXSKIP; If more than ''MAXSKIP" consecutive cores cannot befitted, the iterative fitting process is 
terminated. 

THRESH gives the minimum peak amplitude of cores to be fitted by the GAUSSCLUMPS algorithm. The 
supplied value is multipled by the RMS noise level before being used. 

NPAD: The algorithm will terminate when ''NPAD " consecutive cores have been fitted all of which have 
peak values less than the threshold value specified by the "THRESH " parameter. (From the source code 
cupidGauss Clump s.c, one can see that the algorithm will do the same thing when "NPAD " consecutive 
cores have pixels fewer than "MINPIX".) 

MAXBAD: The maximum fraction of bad pixels which may be included in a core. Cores will be excluded 
if they contain more bad pixels than this value. 

VELORES: The velocity resolution of the instrument, in channels. 
MODELLIM: Model values below ModelLim times the RMS noise are treated as zero. 
MINPIX: The lowest number of pixel which a core can con tain. 
FWHMBEAM: The FWHM of the instrument beam, in pixels. 

MAXCLUMPS: The upper limit of the cores to be fitted. Set to a large number so this parameter do not 
take effect. 

MAXNF: The maximum number of evaluations of the objective function allowed when fitting an individual 
core. Here it is just set to a very large number to guarantee all the cores to be fitted. 
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Table 2: Parameters used in the CLUMPFIND fitting of the 13 CO data cube). 



Parameter 


Value 


DELTAT 


5.0*RMS 


FWHMBEAM 


2.0 


MAXBAD 


0.05 


MINPIX 


16 


NAXIS 


3 


TLOW 


5*RMS 


VELORES 


2.0 



Table 3: Parameters used in the GAUSSCLUMPS fitting of the extinction map. 



Parameter 


Value 


WWIDTH 


2 


WMIN 


0.01 


MAXSKIP 


50 


THRESH 


5 


NPAD 


50 


MAXBAD 


0.05 


VELORES 


2 


MODELLIM 


0.05 


MINPIX 


16 


FWHMBEAM 


2 


MAXCLUMPS 


2147483647 


MAXNF 


200 
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Table 4:: Properties of the cores found in 3D 13 CO 
data cube. 
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Table 4 — continued from previous page 
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Properties of the 10 most massive CO cores. Following the identifiers in column 1, the next two 
columns are the RA and Dec of the center of the fitted 13 CO cores. The semi-major axis -Rmajor and semi- 
minor axis -Rminor of the fitted 13 CO cores follow. The sixth column is the position angle of 13 CO cores 
(angle from north to the major axis). The next column is the peak temperature of the 13 CO cores. The 
mass M and the virial mass M V i r follow. Column 10 gives the full width to half maximum line width of the 
13 CO line. The last column shows the estimated mean H 2 density of the 13 CO cores. 

Table 5:: Properties of the cores found in the 
smoothed 13 CO data cube. 
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Table 5 — continued from previous page 
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The columns are arranged as those of table HJ 
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Table 6:: Properties of the cores found in extinc- 
tion map. 
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Properties of the 10 most massive cores found in 2MASS extinction maps. Following the core identifiers, 
the next two columns are the RA and Dec of the center of the fitted cores, the the semi-major i? ma jor and 
semi-minor i? m i n0 r of the fitted cores. The sixth column is the position angle of cores (angle from north to 

the major axis). The last column is the core mass M. 



